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Abstract We investigate magnetic properties of an ultracold Fermi gas with pop- 
^rj ulation imbalance. In the presence of population imbalance, the strong-coupling 

theory developed by Nozieres and Schmitt-Rink (which is frequently referred to 
as the NSR theory, or Gaussian fluctuation theory) is known to give unphysical re- 
sults in the BCS-BEC crossover region. We point out that this problem comes from 
how to treat pseudogap effects originating from pairing fluctuations and many- 
body corrections to the spin susceptibility. We also clarify how to overcome this 
problem by including higher order fluctuations beyond the ordinary T-matrix the- 
ory. Calculated spin susceptibility based on our extended T-matrix theory agrees 
well with the recent experiment on a 6 Li Fermi gas. 
PACS numbers: 03.75.Hh, 05.30.Fk, 67.85.Lm. 



1 Introduction 

High controllability of ultracold Fermi gases enables us to study various many- 
body phenomena in a systematic manner 1 - 2 . Using a tunable pairing interaction 
associated with a Feshbach resonance, we can study properties of Fermi super- 
fluids from the weak-coupling regime to the strong-coupling limit in a unified 
manner 3 . Introducing population imbalance to a superfluid Fermi gas^, one can 
simulate Fermi superfluids discussed in various research fields, such as metal- 
lic superconductivity under an external magnetic field, and color superconductiv- 
ity in a dense quark matter 6 . The population imbalance naturally induces mis- 
match of Fermi surfaces between (pseudo)spin-^ and (pseudo)spin-J, components, 
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which suppresses the conventional s-wave superfluid stated. In this case, other 
pairing states, such as the Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) state 8 ' 9 and 
the Sarma state, have been predicted—. 

When we deal with both the tunable interaction and population imbalance at 
the same time, pairing fluctuations are crucial in determining the phase diagram in 
terms of the interaction strength, polarization, and temperature. Even in the weak- 
coupling regime, it has been pointed ou t n i 12 that the FFLO phase is unstable 
against pairing fluctuations in the absence of a lattice potential. However, despite 
the importance of strong-coupling effects, polarized Fermi gases have been so far 
mainly examined within the mean-field level, so that roles of pairing fluctuations 
in this system have not been fully understood yet. One of the crucial reasons for 
this current situation is the breakdown of the strong-coupling theory developed by 
Nozieres and Schmitt-Rink (NSR) 13 in the presence of population imbalanc e 14 > 15 . 
While the NSR theory has succeeded in explaining the BCS (Bardeen-Cooper- 
Schrieffer)-BEC (Bose-Einstein condensation) crossover physics in unpolarized 
Fermi gases, it gives negative spin susceptibility in the presence of population 
imbalance 14 , which is, however, forbidden thermodynamically— . 

To resolve the above mentioned problem that the NSR theory possesses, we 
investigate a polarized Fermi gas in the BCS-BEC crossover region. We point out 
that one must carefully treat the pseudogap phenomenon associated with pairing 
fluctuations, as well as many-body corrections to the spin fluctuations in the pres- 
ence of population imbalance. To overcome this problem, the ordinary T-matrix 
theory (which has been also frequently applied to the unpolarized case) is found to 
be still not enough. Further extending the T -matrix theory to include higher order 
fluctuation effects, we obtain the expected positive spin susceptibility in the entire 
BCS-BEC crossover region. The calculated spin susceptibility agrees well with 
the recent experiment on a 6 Li Fermi gas. Using this extended T-matrix theory, 
we also determine the superfluid phase transition temperature T c in a polarized 
Fermi gas. 

2 Formulation 

We consider a two-component Fermi gas with population imbalance (N+ > 
where N is the number of atoms with pseudospin a =t,4-)- Since recent exper- 
iments are all using a broad Feshbach resonance, we employ the single-channel 
BCS model, described by the Hamiltonian, 

H = £ %P,o c p,o c P,o - U £ C p+q/2,t C -p+q/2,i C -P , +9/U C P'+9M- W 

p° p,p',q 

Here, cj, jCT is the creation operator of a Fermi atom with momentum p and pseu- 
dospin a. §p iCT = e p Ho = p 2 /2m — ji a is the kinetic energy, measured from the 
chemical potential jj, a (where m is an atomic mass). When we write the chemical 
potential as fi a = fi + oh, the Hamiltonian in Eq. (fl]) effectively describes a system 
under an external magnetic field h. The pairing interaction —U (< 0) is assumed 
to be tunable by a Feshbach resonance, which is related to the s-wave scattering 
length a s as A%a s jm = -U I [l -UHfim/p 2 )] (where m c is a cutoff energy). In 
this scale, the weak-coupling BCS regime and the strong-coupling BEC regime 
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are, respectively, characterized by (kpa s ) <, — 1 and (kpa s ) > 1 (where 
is the Fermi momentum). The region —1 < (kpa s )~ l < 1 is called the crossover 
region. 

Strong-coupling effects are conveniently described by the self-energy E p (i(O n ) 
in the single-particle thermal Green's function, 

G p ,a[i(On) = [G° pa (i(O n )- 1 -Z p , o (i(0 n )}-\ (2) 

where ft),, is the fermion Matsubara frequency, and G po (i(0„) = [i(Q n — E, P ,o]~ l is 
the non-interacting Green's function. In this paper, we take the self-energy as 

Zp,o(i(D n ) = T £ r q (iV n )G q - p ,- a (iv n - i(O n ). (3) 
qjv„ 

Here, V„ is the boson Matsubara frequency. The particle-particle scattering matrix 
r q (iv n ) = —U /[l — U n q (iv n )] involves fluctuation contributions described by the 
sum of ladder-type diagrams 17 , where 

n g (iv n ) = T £ G p+q/2t (iv n + iO) n )G% +q/2 l (-ico n ) (4) 

p,i(0 n 

is the lowest-order pair-propagator. 

Th self-energy in the ordinary T-matrix theory is given by Eq. (O where 
G p .o{i(O n ) is replaced by G° p(y (i(O n ), as 

£™ A M.) = T £ r q (iV n )G° q _ Pi _ a (iV n - i(On). (5) 
q,iv„ 

The NSR theory also uses Eq. (01, where Eq. (0 is expanded to 0(Z™ A (i'ft)„)), 
as 

G N p f(ico n ) = G%(i(o n ) + Gl (i(O n )Z™ K {m n )Gl c {i(O n ). (6) 

For clarify, we call the strong-coupling theory using the self-energy in Eq. OJ the 
extended T-matrix theory in this paper. 

The superfluid phase transition temperature T c is determined from the Thou- 
less criterion 

r 9 - 1 (/v„ = o) = o. (7) 

When the highest T c is obtained at q = 0, the ordinary uniform superfluid phase 
is realized. The FFLO state, on the other hand, corresponds to the case of q ^ 0. 
However, since the latter state is known to be unstable against pairing fluctuations 
in the absence of a background lattice potentia l 1 ^ 12 , we only consider the case of 
q = in this paper. As usual, we solve Eq. (0), together with the number equations 

N °=U = T L GpA^O- (8) 

p,i(o n 
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Fig. 1 Feynman diagrams describing the spin susceptibility %. (A) ^etma: extended T-matrix 
approximation. (B) ^tma : r -matrix approximation. (C) ^nsr : NSR theory. The double lines 
and single solid lines describe the full Green's function G and the free Green's function G°, 
respectively. We note that, the so-called 'Aslamazov-Larkin' diagrams vanish identically in the 
spin susceptibility. (D) scattering matrix F in the T-matrix approximation, where the dashed 
lines describe the interaction —U. 



3 Spin susceptibility and many-body corrections 



The spin susceptibility % ls given by 



X = 



h 



p,ia>„,< 



dG pa (i(Dn) 



dh 



ft->0 



(9) 



Substituting the Green's function given by Eq. (© into Eq. (0, we automatically 
obtain, not only self-energy corrections, but also vertex corrections to the spin 
susceptibility in a consistent manner. In the extended T-matrix theory (where the 
self-energy in Eq. (O is used), many-body corrections to the spin susceptibility 
(= Xetma) are diagrammatically given by the first line in FigQ] Here, the first term 
(Al) involves many-body corrections to % through the density of states p( ft)). That 
is, noting that the spin susceptibility % is deeply related to p(co) at the Fermi level 
(ft) = 0), when pairing fluctuations induce the pseudogap in p(ft)) above T c , % is 
also suppressed due to the suppression of p(ffl) around ft) = 0. To understand the 
role of (A2), (A3) and higher order diagrams in FigQ] it is convenient to approxi- 
mate r to the lowest-order contribution —U. Then, the sum of these diagrams are 
found to be similar to the random-phase approximation (RPA), as 



Xetma — 



X 



(10) 



where % is given by the diagram (Al) in FigQ] 

When one uses the NSR Green's function G NSR = G° 



G°i:™ A G in Eq. 



(O, the resulting spin susceptibility (= Jnsr) ls diagrammatically described by 
the third line in FigQ] Reflecting that the self-energy is only taken into account 
to 0(E™ A ), the dia gram (Al) is now approximated to (CI) and (C2). In this 
regard, we note that the NSR Green's function G NSR is known to overestimate 
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Fig. 2 (Color online) Spin susceptibility % at T c , as a function of the interaction strength. Xo 
is the spin susceptibility in a free Fermi gas at T = 0. The interaction is measured in terms 
of the inverse scattering length a s . The inset shows Xtma magnified around the region where 
Xtma<0. 



the magnitude of the pseudogap 17 , leading to the unphysical negative density of 
states around the Fermi level in the BCS-BEC crossover region. Because of this, 
the suppression of the spin susceptibility by the pseudogap effect is also overes- 
timated in #nsr, leading to the negative spin susceptibility, as shown in Figf2] In 
the strong-coupling BEC limit, one findsi& 



Xnsr = Xo — ( ^ J Mj)^ 



(11) 

h=0 



where Xo = (27") 'X^cosh 2 (^ p /2T) is the susceptibility of a free Fermi gas. 
N° = Y.p[e^ Ep ~ fl ~ h} + l]" 1 and r c BEC = 0.218r F is the superfluid phase transition 
temperature in the BEC limit (where is the Fermi temperature). The second 
term in Eq. dill) just arises from the diagram (C2) in Fig. [T] which is related to 
fluctuation corrections to the density of states^. As shown in Fig|2j this term is 
the origin of the negative spin susceptibility in the BEC regime. 

The problem of the negative density of states is absent in the ordinary T-matrix 
theory (where Eq. (O is used) 17 . In this case, the (Al)-type diagram in FigQ]is 
correctly included as (Bl) in the second line of Fig|2j However, when we evaluate 
X based on this theory {= #tmaX although the situation becomes much better 
than the NSR result, #tma still becomes negative for 0.5 <5 (kpas) -1 <, 1.1. (See 
the inset in FigO. This is because the RPA-type series in the extended T-matrix 
theory is truncated at the first order ((B2) in FigO in the ordinary T-matrix theory. 
Evaluating #tma by approximating F to — U, one obtains 

Xtmk-X-XoUX- (12) 

While Eq. dlOl ) is always positive, Eq. dl2b becomes negative when UXo > 1- 

The above analyses clearly indicate that one needs to carefully treat the pseu- 
dogap effect and the RPA-type many-body corrections to the spin susceptibility 
in the presence of population imbalance. Since they are both correctly treated in 
our extended T-matrix theory, the calculated spin susceptibility #etma is positive 
in the entire BCS-BEC crossover. As shown in Figf2l #etma gradually decreases 
with increasing the interaction strength, and almost vanishes when (^f^s) 1 ^ 1, 
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Fig. 3 (Color online) Comparison of calculated spin susceptibility Xetma (solid line) with the 
recent experiment on a 6 Li Fermi gas (where the spin susceptibility is measured by using in situ 
imaging of dispersive speckle patterns) (solid circles)^ .. In this experiment, the temperature is 



fixed at the value of T c at (kp a s 
superfluid state. 



= —0.35. Thus, the right side of the vertical line is in the 




Fig. 4 (Color online) Superfluid phase transition temperature T c of a Fermi gas with population 
imbalance. 



because the spin degrees of freedom almost dies out due to the formation of tightly 
binding singlet pairs in the BEC regime. 

We compare our theoretical result with the recent experiment on 6 Li 20 in Figf3] 
The calculated spin susceptibility Xetma well agrees with the observed suscepti- 
bility in the normal state (the left side of the vertical dotted line). We emphasize 
that there is no fitting parameter in obtaining Figf3] 

So far, we have discussed spin susceptibility, which is obtained in the zero 
polarization limit (h — > 0). Here, we briefly consider the case of finite polariza- 
tion. Figure |4] shows the phase diagram of a polarized Fermi gas in terms of the 
interaction strength, polarization P = (N+ — A^)/(A^ +A^), and the temperature. 
In the weak-coupling BCS regime, T c vanishes at a critical polarization P c . The 
magnitude of P c increases, as one passes through the BCS-BEC crossover region. 
In the BEC regime, one finds T c = r c BEC (l -P) 2 / 3 . 

We briefly note that, although we have only considered the second order phase 
transition by using the Thouless criterion in FigJH the first order transition may 
exist below P c 10 i 15 . To examine the first order phase transition of a polarized Fermi 
gas, one needs to take into account the possibility of phase separation between the 
paired molecules and fully polarized Fermi atoms, which remains as our future 
problem. 
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4 Summary 

To summarize, we have investigated an ultracold Fermi gas with population im- 
balance. We clarified that the breakdown of the NSR and T-matrix theories in the 
presence of population imbalance can be eliminated by carefully treating the pseu- 
dogap effects and the RPA-type many-body corrections to the spin susceptibility. 
The resulting extended T-matrix theory gives the expected positive spin suscepti- 
bility in the entire BCS-BEC crossover region. The calculated spin susceptibility 
well agrees with the recent experiment on a 6 Li Fermi gas. Since polarized Fermi 
gases are deeply related to metallic superconductivity under an external magnetic 
field, our results would be useful for the further development of the physics of 
polarized Fermi superfluids. 

Acknowledgements We would like to thank S. Watabe, Y. Endo, D. Inotani, and R. Hanai for 
useful discussions. Y. O. was supported by Grant-in-Aid for Scientific research from MEXT in 
Japan (22540412, 23104723, 23500056). 



References 

1. S. Giorgini, L. Pitaevskii, and S. Stringari, Rev. Mod. Phys. 80, 1215 (2008). 

2. I. Bloch, J. Dalibard, and W. Zwerger, Rev. Mod. Phys. 80, 885 (2008). 

3. C. Chin, R. Grimm, P. Julienne, and E. Tiesinga, Rev. Mod. Phys. 82, 1225 
(2010). 

4. M. W. Zwierlein, A. Schirotzek, C. H. Schunck, and W. Ketterle, Science 311, 
492 (2006). 

5. G. B. Partridge, W. Li, R. I. Kamar, Y.-A. Liao, and R. G. Hulet, Science 311, 
503 (2006). 

6. W. V. Liu and F. Wilczek, Phys. Rev. Lett. 90, 047002 (2003). 

7. A. M. Clogston, Phys. Rev. Lett. 9, 266 (1962). 

8. P. Fulde and R. A. Ferrell, Phys. Rev. 135, A550 (1964). 

9. A. I. Larkin and Y. N. Ovchinnikov, Sov. Phys. JETP 20, 762 (1965). 

10. For review, see D. E. Sheehy and L. Radzihovsky, Ann. Phys. (N.Y.) 322, 1790 
(2007). 

11. H. Shimahara, /. Phys. Soc. Jpn. 67 1872 (1998). 

12. Y. Ohashi, /. Phys. Soc. Jpn. 71, 2625 (2002). 

13. P. Nozieres and S. Schmitt-Rink, /. Low Temp. Phys. 59, 195 (1985). 

14. X.-J. Liu and H. Hu, Europhys. Lett, 75, 364 (2006). 

15. M. M. Parish, F. M. Marchetti, A. Lamacraft, B. D. Simons, Nat. Phys. 3, 124 
(2007). 

16. G. L. Sewell, Quantum Mechanics and its Emergent Macrophysics, Princeton 
University Press, 2002. 

17. S. Tsuchiya, R. Watanabe, and Y. Ohashi, Phys. Rev. A, 80, 033613 (2009). 

18. T. Kashimura, R. Watanabe, and Y. Ohashi, in preparation. 

19. A. A. Varlamov, G. Balestrino, E. Milani, and D. V. Livanov, Adv. Phys. 48, 
655 (1999). 

20. C. Sanner, E. J. Su, A. Keshet, W. Huang, J. Gillen, R. Gommers, and W. 
Ketterle, Phys. Rev. Lett. 106, 010402 (2011) . 



